Quantum vortices in optical lattices 
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A vortex in a superfluid gas inside an optical lattice can behave as a massive particle moving in 
a periodic potential and exhibiting quantum properties. In this Letter we discuss these properties 
and show that the excitation of vortex motions in a two-dimensional lattice can lead to striking 
measurable changes in its dynamic response. It would be possible by means of Bragg spectroscopy 
to carry out the first direct measurement of the effective vortex mass, the pinning to the underlying 
lattice, and the dissipative damping. 
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The understanding of the static and dynamical behav- 
ior of vortices has been crucial to describe numerous dif- 
ferent situations in superfluids ranging from liquid He- 
lium to high-temperature superconductors [j], 0] . These 
defects can be created by means of an applied magnetic 
field in superconductors and by putting the sample into 
rotation in superfluid Helium, or they can be thermally 
excited in low-dimensional systems where the unbinding 
of vortex-antivortex pairs is at the core of the Berezinskii- 
Kosterlitz-Thouless transition. Low-dimensional super- 
conductors and in particular Joscphson Junction Arrays 
(J JAs) have been for many years a natural playground for 
studying classical and quantum properties of vortices [|| . 
A vortex in a JJA behaves as a massive particle moving 
in a periodic potential and subject to dissipation 0], and 
under appropriate conditions vortices can show quantum 
properties such as interference or tunneling. Among the 
most interesting experiments performed with vortices in 
J JAs we mention the observation of ballistic motion [B| , 
the measurement of the Aharonov-Casher effect for a vor- 
tex going around a charge 0, and the Mott- Anderson 
insulator of vortices Q . 

Optical lattices for atomic gases, which currently are 
under intense investigation [H, 0, [13] , can behave as tun- 
neling junction arrays. In this Letter we analyze vortex 
excitations in an optical lattice and show that a super- 
fluid gas in an optical lattice offers a unique opportunity 
for a direct measurement of vortex properties (such as 
the mass, the coupling to its environment or the pinning 
potential) via a Bragg spectroscopy experiment. This is 
in contrast to the JJA case, where only indirect mea- 
surements based on transport properties are available. 
The Bragg spectroscopy technique [U, [H, [l3| has been 
appealed to for a variety of experiments on ultra-cold 
atomic gases, and in optical lattices has been consid- 
ered for a measurement of the excitation spectrum of 
a Bose gas in the Mott-insulator phase 14|, and of its 
coexistence with a superfluid phase in a dishomogeneous 
cloud 0. 

We consider a Bose gas in the superfluid phase inside 
a lattice [Til . 17 1, in a regime where the hopping and 



the on-site repulsion between the bosons are competitive. 
Quantum fluctuations due to the interplay of the local re- 
pulsions and of the hopping have dramatic consequences 
for vortex dynamics. In this case a vortex behaves as 
a macroscopic quantum particle, moving in a periodic 
potential with a mass that we evaluate and show to be 
directly measurable by Bragg spectroscopy. At variance 
from other recent studies of vortices in frustrated optical 
lattices [2, lj| 2(J 21 , 2^, 23| , we discuss the dynamical 
properties of an individual vortex. In order to achieve 
this regime one can either apply a very low frustration 
by means of a rotation of the lattice [2J] generating only 
a few and very weakly interacting vortices, or create a 
vortex excitation by means of phase inprinting [25L l26| . 
Of particular relevance is the very recent observation of 
vortex pinning in co-rotating optical lattices by Tung et 
al. [2J], which indicates that what we propose here is 
within reach of experimental verification. 

The model - We consider a Bose gas at zero temper- 
ature inside a square lattice with lattice constant a and 
N s lattice sites. We assume that the system can be de- 
scribed by a single-band Bose-Hubbard Hamiltonian 27 1 



H.c. 
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where b\ and bi are the creation and annihilation opera- 
tors for a boson on the i-th site and fii = b[b% is the num- 
ber operator. The coupling constant U describes the local 
interaction between bosons, [i is the chemical potential, 
and J the matrix element for hopping between nearest- 
neighbors sites. The on-site interaction energy and the 
hopping energy are given by U — g / G?r |wo(r — R;)| 4 , 
J = -(h 2 /2m)- 1 J driver -R.i)V 2 w (r -fy), in terms 
of the Wannier function w (r) (R; is the coordinate of 
the i-th site). Here g = 4:nh 2 a sc /(\/2Trml±) is the repul- 
sive interaction strength in the case where the transverse 
size l± of the lattice is larger than the scattering length 

^se- 
lf the average number of bosons per site is much larger 
than one, the field operators can be approximated as bi ~ 
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V^exp (i<f>i), with <j>i being the phase operator on the i- 
th site. The Bose-Hubbard model can be recast into the 
quantum phase Hamiltonian 



H = -JnJ2' 



tfj^n? (2) 



where /2 = 2U — fj, — 1. The number operator has been 
expressed in terms of the fluctuations around its average 
value n, hi = n + Shi. The number fluctuation opera- 
tor and the phase are canonically conjugate variables, 
[5h i ,e ±l ' l ' j ] = Sije ±l ^ j . The regime that we consider 
throughout this work is Jh ^> U: the system is deep 
in the superfluid region, but quantum fluctuations are 
present and play a crucial role in the vortex dynamics. 

Vortex properties - The presence of a static vortex 
inside the lattice can be described to a good approxi- 
mation by a phase distribution of the boson field given 
by 



arctan 
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where x, y are the coordinates of the center of the vortex. 
Deep in the superfluid regime and at temperatures much 
lower than Ju/Kb, phase rigidity ensures that again to 
a good approximation, a moving vortex can still be de- 
scribed by Eq. © but with a time-dependent position 
of the vortex center. The existence of a vortex mass can 
be understood qualitatively by noting that if a vortex 
moves of a distance of the order of a in a time St = a/v, 
v being its velocity, the phase difference 5<f>ij at each bond 
changes in time as 6cj>ij — 4>ij (t+a/v) — cf>ij (t). Due to the 
commutation relation between the number and phase op- 
erators, a time-dependent phase leads to a contribution 
to the energy which is quadratic in the vortex velocity 
(see the second term of the r.h.s of Eq.©). The problem 
of calculating the vortex mass can then be reduced to 
find the phase differences across junctions at times t and 
t + a/v. 

An effective action for a vortex in a lattice can then 
be obtained by inserting Eq. (J3]) in the Bose-Hubbard 
model in Eq. j2]) and then expressing the resulting action 
in terms of the vortex coordinates r(i) = {x(t),y(t)} 
The on-site repulsion term provides a kinetic energy term 
T = (M v /2)r 2 , where the vortex mass in a lattice of size 
L is 



(4) 
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The vortex mass thus scales linearly with the boson mass, 
increases with the width w of the Wannier function, and 
decreases with the scattering length. 

Let us consider for an illustration the 87 Rb lattice real- 
ized by Greiner and coworkers [13] , in the case Vq = AE r 
where Vq and E r are the well depth of the optical lat- 
tice and the recoil energy respectively. For such a system 



in 2D we have a — 426 nm, w ~ 96 nm, l± = 5 /xm, 
a sc = 5.5 nm and L = 75 /xm, and the vortex mass is 



M v ~ 29mln(L/a) ~ 150m : 
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The behavior of the vortex mass as a function of the 
potential well depth (which affects w) and of the size of 
the lattice is depicted in Fig. [1] 
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FIG. 1: The vortex mass M v (in units of the boson mass m) 
as a function of Vo/E r and L/a, for the case a — 426 nm, 
l± — 5 fim, and a sc — 5.5 nm. 



The effective potential seen by the vortex has been 
numerically evaluated in the context of JJAs and, in 
the case of a vortex moving in the x direction inside 
a large two-dimensional array the effective potential is 
periodic [28| . 



U v (x) = 0.1Jn[cos(27ra/a) — 1] 



(6) 



In the presence of a vortex the whole array thus behaves 
as a macroscopic particle of mass M v moving in a periodic 
potential. For such a macroscopic object one has to also 
take into account the interaction with the environment, 
and the main source of damping is due to the interaction 
with the long-wavelength phase modes that are excited 
during vortex motion. This damping has been analyzed 
in detail in the context of JJAs (see for example 2^, 30j|). 
In this Letter for simplicity we shall just comment on how 
our results are modified in accounting for dissipation. 

Dynamic structure factor - We turn to a calculation 
of the dynamical response of the Bose gas inside a lat- 
tice, both in the absence and in the presence of a vortex. 
The central quantity of interest is the dynamic structure 
factor, which in a tight-binding scheme takes the form 



S(q, w) = J alte 1 ' 



iq-(Ri-Rj; 



(Shi(t)5h 3 (0)). (7) 



This spectrum can be measured in experiments of Bragg 
spectroscopy ll|, |l2j: two probe laser beams, with fre- 
quencies u)% and LJ2 = u>i + u> and wave-vectors ki 
and k2 = ki + q, scatter on the boson gas and the 
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spectrum measures the probability of momentum trans- 
fer hq at energy hui [l3|. The /-sum rule gives the 
first spectral moment Mi(q) as M\ = J S(q,uj)ui duj = 
7p-(0|[5n q , [H , Jn q ] |0) , <5n q being the Fourier transform 
of 5hi. 

We first consider the case in which no vortex is present. 
Inside the superfluid regime ( Jn 3> U), it is enough to 
consider long- wavelength phase fluctuations, as described 
by expansion of the cosine in the phase Hamiltonian up 
to second order. In this limit the Hamiltonian is eas- 
ily diagonalizcd in Fourier space by means of the trans- 
formations </> k = [UN s /(hfiu)] 1/2 (a k + aL k ) /V2 and 

n k = (N s hVL^/U) 1/2 fa k - a)_^j /iy/2, with the result 

H = J2 m * («k«k + (8) 

where Q 2 , = (2JnU/h 2 )[2—cos(k x a) — cos(k y a)]. Here the 
quasi-momentum k = [k x , k y ) is inside the first Brillouin 
zone. By taking into account the time dependence of 
the particle number fluctuation operator dictated by Eq. 
© , it is straightforward to obtain the dynamic structure 
factor as 

<?(q )W ) = Mks(u;-Q q ). (9) 

The physical interpretation of Eq. is clear: small-g 
absorption occurs at a frequency corresponding to the 
dispersion relation of the Goldstone sound mode in the 
lattice. In this case Mi(q) = (Jn/h)[2 — cos(q x a) ~ 
cos(g„o)]. 

The presence of a vortex can induce, besides changes 
in the sound wave spectrum, specific contributions asso- 
ciated with excitations of vortex motions. Several differ- 
ent situations can be envisaged for the latter, but here 
we discuss the interesting case in which the hopping pa- 
rameter J is sufficiently large that the vortex is pinned 
to a minimum of the periodic potential given by Eq. 
In this case the vortex dynamics is associated with small 
oscillations around its equilibrium position and can be 
described by the harmonic oscillator Hamiltonian 

H v = hl v r 2 + hl v n 2 v r 2 , (10) 

where we have defined Q„ = (0.1Jn/M v ) 1 / 2 2n/a. By 
noting that 

U 2 ■ * 

<*fti(ti)toj(*a)) = Tj(^(tiWj(tj)) , (ii) 
h 

using the expression given in Eq. for the phase dis- 
tribution and performing an average over the vortex de- 
grees of freedom with the help of Eq. (flQ|) , it is possible 
to write the contribution of the vortex to the dynamic 
structure factor as 

S v (q lU ) = ^ 2 6(u - n v ) . (12) 



Instead of a q-dependent resonance as in Eq.®, the 
Bragg scattering acquires a resonance at a well defined 
frequency f2„ indicating that the whole lattice responds 
collectively having the properties of a single macroscopic 
particle, the vortex. This is the main result of this pa- 
per. Under the conditions specified above, the presence 
of a vortex induces a resonance at a frequency that al- 
lows access to the vortex mass. Let us remark that this 
resonant behavior is related to the presence of the lattice 
and to the existence of quantum fluctuations originating 
from the local repulsion. The Bragg spectrum of a vortex 
in a Bose-Einstein condensate is otherwise determinated 
by a dispersion relation [l3j |. The peculiar dependence 
of the spectral strength in Eq. (TT2")) on the transferred 
momentum q is due to the coupling between the exciting 
radiation and the lattice: at low momentum all phases 
in the lattice are excited and the dynamic response is 
enhanced. We finally should comment on the fact that 
Eq. (fT2|) does not fulfill the /-sum rule: this should come 
as no surprise, as this expression is valid only at low en- 
ergy. 

The coupling to long-wavelength phase fluctuations 
provides the main dissipation mechanism for the vortex 
motions 29,30,31]. To a first approximation this results 
in Ohmic damping on the vortex. In the presence of dis- 
sipation the delta function in the dynamical response is 
smeared and acquires a finite width proportional to the 
dissipation strength. 

In summary, in this Letter we have discussed some 
main aspects of quantum dynamics of a vortex in an 
atomic superfluid gas inside an optical lattice. We have 
specifically considered the situation in which the vor- 
tex is pinned by the lattice potential and only executes 
small oscillations around its energy minimum. In this 
case Bragg spectroscopy should allow a measurement of 
the vortex mass. One can envisage other situations in 
which to study vortex dynamics. Experiments on quan- 
tum tunneling/coherence of vortices seem to be within 
experimental reach. 
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